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ABSTRACT 

Stellar photometry derived from the INT/WFC Photometric Ha Survey of the North- 
ern Galactic Plane (IPHAS) can be used to identify large, reliable samples of A0-A5 
dwarfs. For every A star, so identified, it is also possible to derive individual reddening 
and distance estimates, under the assumption that most selected objects are on or near 
the main sequence, at a mean absolute r' magnitude of 1.5 - 1.6. This study presents 
the method for obtaining such samples and shows that the known reddenings and 
distances to the open clusters NGC 7510 and NGC 7790 are successfully recovered. A 
sample of over 1000 A stars is then obtained from IPHAS data in the magnitude range 
13.5 < r' < 20 from the region of sky including the massive northern OB association 
Cyg 0B2. Analysis of these data reveals a concentration of ^200 A stars over an area 
about a degree across, offset mainly to the south of the known 1-3 Myr old OB stars 
in Cyg 0B2: their dereddened r' magnitudes fall in the range 11.8 to 12.5. These are 
consistent with a ^ 7 Myr old stellar population at DM = 10.8, or with an age of 
~ 5 Myr at DM = 11.2. The number of A stars found in this clustering alone is 
consistent with a lower limit to the cluster mass of ^ 10'' Mr. 
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ciations: individual (Cyg 0B2, NGC 7510, NGC 7790) 



open clusters and asso- 



1 INTRODUCTION 

In the spectra of stars, the Ha line is a prominent marker 
both as an emission line and as an absorption feature. Ha 
surveys are generally seen eis vehicles for identifying emission 
line nebulae and stars associated with, primarily, the early 
and late phases of stellar evolution. In the case of IPHAS, the 
INT/WFC Photometric Ha Survey of the Northern Galac- 
tic Plane (Drew et al. 2005), the quality of the digital pho- 
tometry in the 3 survey filters, r', i' and Ha, also allows 
the strong Ha absorption associated with A-type stars near 
or on the main sequence to be used as a tool to select these 
early-type objects easily and reliably. Furthermore, the pho- 
tometry is sufficient to provide estimates of reddenings and 
distances. The aim of this paper is to demonstrate how this 
capability may be exploited. In so doing, we shall verify the 
method against two test cases, and present an example of 
this selection at work. 

We shall show that the selection of A stars from the 
IPHAS (r' — Ha, r' — i') colour-colour plane yields mainly 
samples of A0-A5 dwarf candidates. These early-type dwarfs 



of intermediate mass are expected to remain near the main 
sequence for not much longer than ~100 million years, and 
so will typically trace relatively young environments - which 
are a defining feature of the thin disk population of our 
Galaxy. Objects of this type selected via their maximal Ha 
absorption have been exploited before as probes of structure 
within the galactic plane: in a series of papers, S. McCuskey 
and collaborators identified them in objective prism data 
down to ~ 12 and then mapped their distribution (see 
e.g. McCuskey & Houk 1971, McCuskey & Lee 1976). This 
technique can be picked up again and applied to the IPHAS 
database starting from a bright limit of ~ 13th magnitude, 
reaching down to '--^20th. As IPHAS approaches completion, 
and its catalogue database is becoming public (Gonzalez- 
Solares et al. 2008) , the targeted selection of A stars can be 
applied to a range of problems across the northern Galactic 
plane, ranging from focused studies of individual clusters to 
tracing large scale Galactic disk structure. 

The main application considered here is of the more 
localised type - to the still enigmatic Cyg 0B2 associa- 
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tion; data from 21 IPHAS fields covering tlie slcy in and 
around tliis massive OB association are presented and anal- 
ysed. Cyg OB2 is the most impressive of the northern hemi- 
sphere's Galactic OB associations. Our aim is to use the 
A-star distribution to better define the extent and nature of 
this massive cluster. This example also serves to show that 
even very young regions - 1-3 Myrs old in this case (Hanson 
2003) - can be usefully explored by this means. 

Cyg 0B2 (VI Cygni) is home to the only known north- 
ern hemisphere 03 stars (Walborn ct al. 2002). It was first 
recognised for the very massive cluster that it is by Red- 
dish, Lawrence & Pratt (1967), who carried out the first 
extensive study of the region. They placed it at a distance 
of 2.1 kpc, based on the apparent magnitude of late B stars. 
They also defined its extent as a 29 pc x 17 pc ellipse cen- 
tred on RA 20 32 30 Dec 4-41 30 00 (J2000), and placed its 
mass in the range 0.6 - 2.7x10* Mq. However now, mainly 
as the result of the work by Massey & Thompson (1991), 
the main concentration of OB stars is recognised as lying 
just to the NE of the VI Cyg No 8 trapezium located at 
RA 20 33 06 Dec +41 22 30 (J2000), at the nearer distance 
of 1.7 kpc {DM = 11.2). This was not so clearly evident, 
to begin with, because of the substantial reddening of the 
main association: Ay ~ 6 is typical for Massey & Thomp- 
sons's (1991) OB population. On the basis of 2MASS NIR 
photometry of the region, Knodlseder (2000) argued that 
Cyg OB2 could be even more massive than Reddish et al's 
(1967) estimate, perhaps reaching as much as 10^ M©, and 
could contain as many as 2600±400 OB stars. This has since 
been challenged by Hanson (2003) , who has also shown that 
recent recalibrations of OB star absolute magnitudes may 
require a further reduction of the distance to Cyg 0B2 to 
~1.4 kpc (or DM = 10.8). 

The paper is organised as follows. The next section 
presents the main features of IPHAS and the resultant 
database. Then in Section 3, the IPHAS (r' - Ha,r' - i') 
colour-colour plane is discussed from the perspective of the 
reliable selection of A stars from it, how they may be dered- 
dened, and how this selection maps onto a mean absolute 
magnitude. We then present two test cases, in the form of 
the well-studied open clusters NGC 7510 and NGC 7790, to 
show that the selection works as expected in recovering their 
known distances and reddenings. This prepares the way for 
the challenge of exploring the Cyg 0B2 region in section 5, 
in which the full impact of the early A star absolute mag- 
nitude sensitivity to ago, at the youngest ages, is scon. Wc 
conclude with a comment on our findings and on further de- 
velopment and application of the A-star selection technique. 



2 OBSERVATIONS: IPHAS PHOTOMETRY 

A full description of IPHAS data-taking and data^processing 
has been given by Drew et al. (2005). Here, we mention again 
only the key features of the observations. 

All IPHAS imaging is obtained with the Wide Field 
Camera (WFC), mounted on the 2.5-metre Isaac Newton 
Telescope, in La Palma. The WFC is an imager compris- 
ing 4 AR-coated, thinned 4K x 2K EEV CCDs arranged in 
an L shape, capturing data from an on-sky area of approx- 
imately 0.3 of a square degree. With a pixel dimension of 
13.5 //m, corresponding on-sky to 0.333x0.333 arcsec^, the 
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Figure 1. The positions of the extracted IPHAS field-pair over- 
laps in and around CygOB2. The black concentric circles in be- 
tween fields 5985 and 6010 mark the often-quoted centre of Cyg 
OB2 (e.g. as in Knodlseder 2000). 

instrument is able to exploit the sub-arcsecond seeing fre- 
quently encountered at the telescope site. All the observa- 
tions of the Cyg 0B2 region, presented here, were obtained 
during the clear moonless nights of 8 and 9 August 2004. 
The data on the open clusters NGC 7510 and 7790 were ob- 
tained on grey nights, in stable conditions, respectively on 
26 July and 22 October 2005. 

A complete IPHAS observation of any sky position con- 
sists of two pointings, the second at an ofi'set of 5 arcmin W 
and 5 arcmin S with respect to the first. This strategy both 
enhances the overall photometric quality of the final survey 
and ensures that stars falling into a gap between the mo- 
saiced CCDs in one exposure are captured in the partner 
exposure. In this study we only consider point sources that 
are well-detected in both exposures: this limits the sky ex- 
tent of the overlap between each field pair to a near-square 
of side 31 arcmin. Within this there is an irregular grid of 
inter-CCD gaps, with the result that the sample of point 
sources discussed here drawn from each field pair comes from 
a sky area of 0.22 deg'^: since there is very little overlap be- 
tween adjacent field pairs when they are combined in this 
way, all the Cyg 0B2 fields together span close to 90 per- 
cent of the total available area of ~ 5 deg^. The rough sky 
outlines of the IPHAS field-pair overlaps in the Cyg 0B2 
region arc shown in Fig. 1. The coordinates of the centres 
of all catalogue extractions discussed in this paper are listed 
in Table 1, along with other relevant data. 

Throughout this paper field pairs are named just by the 
4-digit number (column 1 in Table 1) that tags the initial 
field pointing in the IPHAS database. 

At each WFC pointing, imaging is carried out through 
three filters. These filters are narrow-band Ha (FWHM 
95 A), Sloan r' and Sloan i' . The exposure times in these 
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Table 1. Co-ordinates of field-pair centres, extracted sky areas and object counts. 



IPHAS Extraction Centre Coordinates Extracted area Star counts: 13.5 < r' < r' 
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(Cyg OB2 region totals:) 910 1250 



three filters are 120 sec (Ha), ,30 sec (r') and 10 sec {i'). 
Since it is required that every star considered is imaged in 
both overlapping pointings, the effective exposure times are 
twice as long. Every night of IPHAS observing incorporates 
observations of standard fields, in order to permit photomet- 
ric calibration in the three observed bands. Presently the 
Ha narrow-band observations are tied to the r'-band cali- 
bration in the same way as described by Gonzalez-Solares 
et al. (2008). 

A further criterion applied in our selection of candidate 
A stars is that the CASU pipeline processing and cataloging 
the data, identifies them as point-like in r' and i' , to a high 
confidence level (i.e. as members of morphology classes - 
1 and -2, Irwin 1985 and 1997). The photometric error in 
the resultant catalogue of point sources ranges from entirely 
negligible at the adopted bright limit (13th magnitude in 
r'), to, typically, 0.03 in r', 0.06 in (r' - Ha) and 0.04 in 
(r' — i') at the chosen faint limit (r' = 20). 



3 METHOD 

3.1 The IPHAS colour-colour plane, illustrated 
by an example from Cyg OB2 

An example of derived colours for a final selection of point 
sources is shown in Fig. 2. These are the data for IPHAS field 



6010, which includes the nominal centre of Cyg OB2. All 
2943 stars plotted fall in the magnitude range 13 < r' < 20. 
To understand what is implied by this source distribution we 
make use of the solar-motallicity synthetic tracks that were 
given in Drew ct al. (2005). On Fig. 2 we have superimposed 
a curve that defines an approximate upper bound on the 
colour-colour plane location of early A dwarf stars. More 
precisely this is the reddening line derived using as template 
the solar metallicity AOV fiux-calibrated spectrum from the 
library due to Pickles (1998), shifted to higher r' — Ha by 
0.03. The reddening law used here is a mean Galactic law, 
with i? = 3.1 (Howarth 1983). 

The main features of the distribution of stars in the 
colour-colour plane (Fig. 2) are (i) a foreground unreddened 
main sequence, begiiiiiiiig at a spectral type of ^GO (at 
r' = 13); (ii) an intermediate group of lightly reddened stars 
{Ay < 3), for which the earliest detected spectral type shifts 
to ~F0 (at r' ~ 14.5); (iii) a concentration of stars spread 
between the reddened Ay = 4.5 and Av = 7 MS lines, and 
reaching down to and below the early-A dwarf limit line. It 
can be deduced from the observed magnitudes and colours 
that the stars in the first two groups will be located at dis- 
tances of no more than ^500 pc. This impression is consis- 
tent with a study by McCuskey and Houk (1971) of A stars 
brighter than y = 13 in which they remarked on a declining 
stellar density between 500 pc and 1 kpc in this part of the 



4 Janet E. Drew et al. 




0.4 0.8 1.2 1.6 2.0 2.4 2.8 
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Figure 2. The {r' — Ha, r' — i') colour-colour plane of IPHAS field 6010 point sources, selected from the magnitude range 13 < r' < 20. 
The black line running along below most of the data points and labelled 'Early A dwarf limit' is the reddening line below which the 
great majority of point sources are expected to be near mains sequence A0-A5 stars. Two reddened main sequences constructed from 
synthetic photometry based on the Pickles (1998) library of flux-calibrated spectra are also drawn (in black also). The red crosses on the 
Ay = 4.5 sequence pick out the data points that make it up (spanning spectral types 05V to M4V). The IPHAS photometric data are 
colour-coded according to source r' magnitude; 13 < r' < 19, darker blue; 19 < r' < 20, light blue. 



Plane, beyond a high concentration found between 200 and 
300 pc (see Fig. 7 of their paper). Group (iii) in Fig. 2 is 
representative of the Cyg OB2 region. Previous work on Cyg 
OB2 itself have claimed typical reddenings in the range from 
4v =i 4 to ~ 7 (Massey &Thompson 1991, Hanson 2003). 
In Fig. 2 it can be seen that the colour-colour plane is well 
occupied in this range. It is worthy of note that redward of 
the Ay = 7 main sequence there is a distinct thinning out of 
stars. This is mainly attributable to the faint limit imposed 
on the included data (cf. Fig. 11 discussed in section 5). 



3.2 The ingredients of IPHAS early-A star 
selection 

At the heart of the method of selection is the choice of a 
suitable cut along the lower edge of the main stellar locus in 
the (r' — Ha,r' — i') plane that separates off candidate early- 
A dwarfs from what will most commonly be later-type stars. 
Once a group of candidate objects has been picked out in 
this way, it is straightforward to use the measured r' — i' of 
each object to derive its reddening, and then to deredden its 
measured apparent r' magnitude. If a suitable averaged ab- 
solute magnitude is available, it becomes possible as a next 
step to determine the distance modulus {DM = 5 log D — 5) 
for the sample. This enables the construction of a redden- 
ing map for the area of sky studied and gives a picture of 



how the A stars are distributed with distance. From this 
approach, one does not expect to derive accurate redden- 
ings and distances for individual candidate A dwarfs, for a 
number of reasons including that the photometry cannot pin 
down the A sub-type to within better than 2 to 3 subtypes, 
and that the absolute magnitude spread at each sub-type 
is in any case significant (Jaschek & Gomez 1998). What 
should emerge is a reliable picture in a statistical sense, in- 
sofar as the mean properties of the selected population of 
objects are well-defined. 

The shape of the lower edge of the main stellar locus in 
the (r' — Ha,r' — i') plane will always follow a segment of 
the reddening curve of an early-A dwarf, as was first noted 
by Drew et al. (2005). Such a curve is therefore appropri- 
ate to use as the cut line. Furthermore, in Fig. 2, it can 
be seen from the labelling of the Ay = 4.5 main sequence 
that spectral types AO to A5 form a line segment running 
essentially parallel to the reddening line labelled "A dwarf 
limit" (actually the AOV reddening line shifted upward by 
A(r' — Ha) = 0.03, in this instance). This is a direct conse- 
quence of the fact that the Ha absorption equivalent width 
reaches a crudely defined maximum of 11 - 13 A, across this 
dwarf sub-type range (cf. comment on this in Drew et al. 
2005, section 5). 

This is not just a quirk of the stars selected by Pickles 
(1998) that were used as a basis for the synthetic photom- 
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etry. To confirm this crucial point, tlie Ha line profiles of 
near-solar metallicity A dwarfs extracted from the Indo-US 
spectral library (Valdes et al. 2004) have been compared 
with those of stars of neighbouring spectral type and lumi- 
nosity class. The results of this are summed up in simplified 
form in Fig. 3. The diagram shows that in the mean, AO- 
5 dwarfs (data on 28 stars combined) are better separated 
from B9 dwarfs (17 stars), than from A6-A8 dwarfs (7 stars). 
Given how spectral classification first arose, this is to be ex- 
pected. The luminosity class comparison is not on such a 
strong statistical footing, as only 16 giant-star spectra make 
up the comparison set. Nevertheless, it is clear that the Ha 
absorption equivalent width in an early-A giant is less, on 
average, than in an early-A dwarf. On an object by object 
basis there is enough scatter that the distinction between, 
especially, >A6 and <A5 dwarfs is not always in the ex- 
pected sense. At the same time there is some compensation 
through the empirical fact that the late A near main se- 
quence range is not as well populated with objects as at 
early A: the bright sample of Abt & Morrell (1995 Table 5) 
shows that around 1 in 3 luminosity class V and IV stars 
are A5 or later. 

This pattern of behaviour translates into reddening 
lines in the (r' — Ha, r' — i') plane as shown in Fig. 4. For 
clarity each reddening line drawn has had the AOV redden- 
ing line subtracted from it, causing the ordinate in the dia- 
gram to become A(r' — Ha). The A0-A5 V reddening lines, 
drawn in black, lie below those for adjacent dwarf spectral 
types and luminosity class HI. The appearance of this di- 
agram gives an indication of the A(r' — Ha) range that 
can be associated exclusively with the early-A dwarfs: it is 
likely to be ~0.02, with the AOV curve used here defining 
A(r' — Ha) — 0.0. Photometric error, together with the in- 
evitable gradual mixing in of (mainly) later type dwarfs as 
the cut line is raised, stands in the way of a fixed recom- 
mendation for placing of the cut line. If the aim of selec- 
tion is to be inclusive of all A dwarfs, choosing a cut at 
A(r' — Ha) — 0.04 or 0.05 would make more sense. Hales 
et al (2008) have obtained spectra of candidates selected us- 
ing this more generous criterion and find that they are the 
expected A stars. 

Here we explore the effect of gradually raising the cut 
on the selection process, to see how increased sample size 
trades against increased likelihood of sample contamination. 
The fit to the unshifted AOV line in the (r' - Ha,r' - i') 
plane that we use is a simple quadratic: 

(r - Ha) = -0.009 + 0.330(r' - i') - 0.0455(r' - i')^ (1) 

The other quantitative issues to settle are the me- 
dian intrinsic r' — i' colour to associate with A0-A5 dwarfs 
(needed for dereddening) , and the median absolute mag- 
nitude, M{r') that should be adopted for them. To asso- 
ciate absolute visual magnitude with dwarf spectral type 
and colour, we use the data based on Hipparcos measure- 
ments presented by Houk et al. (1997). To find suitable me- 
dian values, it is necessary to consider how a typical dwarf 
luminosity function would map onto the A0-A5 range and 
pick out the expected spectral type and colour of the me- 
dian absolute magnitude for this spectral type range. To do 
this, we use the main sequence luminosity function derived 
by Murray et al. (1997) also from Hipparcos data. Houk et 
al. (1997) give Mv = 0.82 for AO, and 1.76 for A5 and note 
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Figure 3. Ha line profiles for spectral types of early-A dwarf 
and adjacent spectral classes. Each profile is the merger of many 
taken from the Indo-US library of coude spectra (Valdes et al. 
2004). The top panel provides a comparison of AO-5 dwarfs (in 
black) with B9V (blue) and A6-8 (red) dwarfs, while the bottom 
shows the difference between AO-5 dwarfs and giants (green) at 
~A type. Before merging, the spectra were rescaled to yield unit 
fiux at 6640 A, to facilitate the comparison without introducing 
further corrections that might alter the very subtle depth differ- 
ences across the central ~80 A of the Ho profile. 



a spread of 0.39 mags in absolute magnitude at each spec- 
tral sub-type. Based on the data in Murray et al. (1997), 
a median My of ^^1.5 is reasonable, for which the intrinsic 
B — V is r^O.l (Houk et al. 1997). From the main sequence 
colours tabulated by Kenyon & Hartmann (1995), intrinsic 
V — Rc — 0.03, and Rc — Ic — 0.06 are implied (spectral 
type A3-4). 

This line of reasoning leads to a median intrinsic r' — i' 
of 0.06, for which the uncertainty is unlikely to be more 
than ±0.02. To derive the colour excess due to extinction, 
E{B — V), from an observed value of r' — i' , this intrinsic 
colour is subtracted to give E{r' — i'), which is then multi- 
plied by 1.58. To obtain the visual extinction. Ay, E{r' — i') 
should be multiplied by 4.91. These transformations apply 
to the case oi R = Av/E{B — V) =3.1, defining the most 
commonly encountered Galactic extinction law. Where non- 
standard reddening laws are thought likely to apply (e.g. R 
approaching 5 as in sightlines dominated by molecular gas) , 
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Figure 4. Reddening lines derived from the IPHAS colour-colour 
plane for spectral types in the early-A dwarf range plotted as 
differences with respect to the AOV reddening line. The same 
colour code is in use as in Fig. 3, e.g. AO-5 dwarf lines are drawn 
in black. The fluxed spectra used to synthesise these tracks were 
taken from the Pickles (1998) library. The lowest line, along the 
horizontal axis, is the AOV reddening line. Note that A{r' —Ha) = 
0.02 runs just below the AO III line, the lowest green line; whereas 
A(r' — Ha) = 0.03, drawn in as a dash-dotted line, falls just 
under the A7V line (in red). The blue line is derived from the 
B9V reddening line. 



it would make more sense to refer all reddenings to either 
the r' or i' band, rather than to V in view of the fact that 
atypical laws vary away from the norm most strongly in the 
blue-to-visual part of the optical spectrum (Cardelli, Clay- 
ton & Mathis 1989). Specifically, the following relations can 
be applied: A^' = 4.13S(r' - i'), and Ai, = 3.13£(r' - i') 
(using tabulated data from Schlegel, Finkbeiner & Davis 
1998). In the cases examined here (sections 4 and 5), previ- 
ous studies have found the average Galactic extinction law 
or laws very like it to be satisfactory. 

The choice of representative absolute magnitude can be 
subject to systematic variation, given that the age spread of 
A dwarfs within 80 pc of the Sun (the samples used by both 
Houk et al. 1997, and Murray et al. 1997) is not guaranteed 
to be the same as that in another potential target region. 
Domingo & Figueras (1999, Fig. 7) have analysed A3-9 main 
sequence stars drawn from essentially the same Hipparcos 
sample as that discussed by Houk et al (1997) and show that 
the earliest among them are indeed quite widely scattered in 
Mv, with mean magnitude ~1.5. But they also demonstrate 
these field stars lie on and smeared above the zero-age main 
sequence (ZAMS), expected to fall at My ~ 2 - that is, these 
stars may have ages of up to ~0.5 Gyr (based on comparison 
with isochrones due to Schaller et al. 1992) . At ages greater 
than ~1 Gyr all stars on or near the early A main sequence 
will have evolved away. 

In Cyg OB2, a region that may be as young as 1-3 mil- 
lion years old (Hanson 2003), we can anticipate that stars 
with the intrinsic (r' — i') colours of early-A dwarfs are more 
massive stars still approaching the ZAMS. Such objects are 
evolving almost horizontally on the HR diagram and are 
likely to arrive on the main sequence as mid-late B stars 
at an age of ~5 Myrs (Siess, Dufour & Forestini 2000). Ac- 
cordingly these are stars with significantly brighter M(r') 
than those destined for the early-A ZAMS, after a pre main 



Table 2. The evolution of M(r') for stars of early-A (r' — «') with 
cluster age. These rough figures are derived from the Z = 0.02 
theoretical isochrones of Siess et al. (2000). Note that at the 
youngest ages of ~3 Myr or less, there is a greatly increased like- 
lihood that these objects will exhibit circumstellar Ha emission, 
rather than maximal absorption (i.e. they will not lie along the 
bottom of the IPHAS main stellar locus). 



age (Myr) 


typical M(r') 


2 


0.0 


3 


0.5 


5 


1.0 


7 


1.4 


10 


1.8 


20-100 


1.9-2.0 


300 


1.6 



sequence evolution taking around 10 Myrs or so (Palla & 
Stabler 1993; Siess et al. 2000). 

Because this age effect can be significant, we give in 
Table 2 some representative values for M(r') as a function 
of age that have been derived from the isochrones of Siess 
et al. (2000). Up to ~10 Myrs, M(r') declines markedly, 
and then remains close to ~2 for a further '^lOO Myrs or 
so - at the end of which, evolution away from the main se- 
quence begins. In the absence of any independent constraint 
on population age or mix, it seems likely that the best guess 
for the general field in the Galactic disk will be around or 
slightly fainter than the Hipparcos value of My — 1.5, or 
M(r') of 1.5 - 1.6. Since the age-dependence has obvious 
and serious ramifications for using the selection of early-A 
dwarfs to analyse Galactic disk sightlines, we begin with two 
control populations to see how well the method works. Our 
guinea pigs are the two open clusters, NGC 7510 and NGC 
7790, aged 6-10 Myrs (Piskunov et al. 2004) and ~100 Myrs 
(Bragalia & Tosi 2006), respectively. 

Finally, it is worth recalling that binarity will also have 
some effect on the average absolute magnitude of a popu- 
lation of early A stars. We make no explicit correction for 
this here beyond favouring M(r') = 1.5, since we estimate 
it amounts to a brightening of the mean in the region of 0.1 
magnitudes. The following considerations suggest this fig- 
ure: roughly a half of all stars are binary systems; the mass 
ratio distribution is not far from fiat (Halbwachs et al 2003); 
hence only in 20-30 % of all binaries is the combined mag- 
nitude likely to brighten by more than 0.5 (see figure 1 in 
Hurley & Tout 1998). 



4 TEST CASES 

4.1 Specifics of the extraction of candidate 
early-A dwarfs 

The steps in the extraction are very simple. First, r' — Ha (as 
a function of r' — i') is used to select stars as candidate early- 
A dwarfs. Once selected, the measured r' — i' is converted 
to a reddening, which then allows the measured r' to be 
corrected to a dereddened magnitude. 

There is no redundancy in our selection of candidate 
objects from IPHAS: we use 3 independent quantities to 
determine 3 properties. Hence it is helpful to bring in e.g. 
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Figure 5. The 6x6 arcmin^ field centred on the test cluster 
NGC 7510, as extracted from the IPHAS database. This is the 
30-sec r' exposures from pointing 7366o. N is up and E to the 
left. 

near-infrared photometry to provide an independent check 
on the selection. For this, we turn to 2MASS J magnitudes - 
where available. Hence, as we select objects below an early- 
A dwarf limit line, we also collect the 2MASS J magnitude 
along with the IPHAS magnitudes and check that the dered- 
dened r' — J colour falls in the early A-dwarf range (set gen- 
erously to be —0.1 < r' — J < 0.3 for A0-5V, cf. Kenyon 
& Hartmann 1995). If the dereddened r' — J for any object 
falls outside this range, the object is dropped from consid- 
eration. This weeding eliminates between a few percent to 
30 percent, depending on the field. If no 2MASS J magni- 
tude is available or it is an upper limit only, we retain the 
candidate A-dwarf unverified, rather than throw it away. If 
we were not to do this, our sample size would be noticeably 
reduced (by a further factor of up to 20-30% in the worst 
cases, again depending on the field - it is greater where the 
reddening is lower). 

In all cases, the bright limit is set at r' = 13.5 in order to 
avoid objects that might be saturated in i' - particularly in 
the reddened fields found in the Cyg OB2 region. The faint 
limit is set at a value appropriate to the sensitivity limit 
for the IPHAS field under consideration. Whilst exploring 
this technique of A-dwarf selection, we treat the position of 
the cut line in the colour-colour diagram as a variable. In 
particular, we obtain samples of candidates for 4 positions 
of the cut line: for A{r'-Ha) = 0.015, 0.02, 0.025 and 0.03, 
where A(r' — Ha) is the upward dispacement of the AOV 
reddening line given in equation 1. 

4.2 The open cluster NGC 7510 

This young open cluster in the northern Galactic Plane has 
recently been studied by Piskunov et al. (2004), along with 
a number of southern examples. The aim of this work was 
to examine the effect of relative youth (main sequence turn- 
on) on the cluster luminosity function. Piskunov et al. con- 
cluded that NGC 7510 is likely to be about 6 million years 
old. Other authors have favoured ages of ~10 million years 
or more (Sagar & Griffiths 1991, Barbon & Hassan 1996, 
Kharchenko et al 2005). The CCD photometric study of the 
cluster by Sagar & Griffiths (1991) earlier found the distance 
modulus to be 12.5±0.3 and that E{B — V) ranges from 1.0 
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Figure 6. The (r' — //o, r' — i') colour-colour plane of the 6x6 
arcmin^ field centred on the young open cluster NGC 7510 (in 
IPHAS fields 7366/o). The point sources plotted are selected from 
the magnitude range 13.5 < i"' < 19.0. The parallel solid and 
dashed lines are the early-A dwarf cut line shifted upward from 
the AOV reddening line by, respectively, 0.03 and 0.015. The piece- 
wise connected line is the E(B — 'V) = 1.2 main sequence con- 
structed from synthetic photometry based on the Pickles (1998) 
library of fiux-calibrated spectra. 



to 1.3, in good agreement with e.g. Lynga's (1988) cluster 
catalogue. 

Our selection of candidate A dwarfs is taken from a 6 x 6 
arcmin'^ region of sky centred on the cluster (Fig. 5). The 
faint limit to the selection is set at r' = 19 - two magni- 
tudes above the pipeline estimate of the magnitude limit for 
the observations of IPHAS fields 7366/o in which NGC 7510 
is located. The lowest and highest cut lines applied to the 
data are shown in Fig. 6, superposed on the colour-colour di- 
agram for NGC 7510 and environs obtained from the IPHAS 
database. 

The results of the selection are shown in Fig. 7. In 
the upper panel, the raw histogram of star numbers per 
AT) Ml. 5 = 0.25 bin is presented. The distance modulus used 
here is derived applying a cluster age-appropriate value for 
M(r') of 1.5 (see Table 2). The extracted sample sizes are 
31, 34, 43 and 53 objects, in order of increasing A(r' — Hd). 
In the middle panel, an arbitrarily-scaled quantity directly 
proportional to the space density of candidate A-dwarfs, pA, 
is plotted againt DM1,5. To arrive at pA, the number of can- 
didate A-dwarfs per distance modulus bin is divided by the 
cube of the bin mid distance: this corrects for the growing 
volume associated with each successive histogram bin. The 
bottom panel shows the relationship between DA/i.s and 
reddening for two samplings of the data. 

The main cluster is very clear, in all three panels of 
Fig. 7. In the histogram (upper panel), only 5 out of the sam- 
ple of 34 stars for A(r' — Ha) — 0.02 (blue line and points) 
fall outside the main peak, centred on DM\.z, — 12.53 ±0.04 
(measured by fitting a single Gaussian to the peak). For 
the largest sample, obtained with A(r' — Ha) = 0.03, 12 
stars out of 53 stray away from the cluster peak (green line 
and points), that itself has shifted downwards slightly to 
DM1.5 = 12.48 ± 0.02. The increased proportion of outliers 
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Figure 7. The tost case of the open cluster NGC 7510. The upper 
panel shows the distribution of derived distance moduli (DM1.5 = 
dereddened r' magnitudes less an adopted median absolute mag- 
nitude of 1.5). The four curves shown are for four choices of posi- 
tioning of the early-A dwarf limit line: A(r' - Ha) = 0.015, 0.02, 
0.025 and 0.03 are respectively in black, blue, red and green. Note 
that as the limit line is raised, in A(r' — Ha) = 0.005 steps, the 
distribution is better populated but also begins to shift slightly 
toward lower mean distance modulus. The middle panel is a con- 
version of the upper panel to an A-star space density distribution, 
obtained by taking the number of stars in each histogram bin plot, 
and dividing by cube of the distance. The lower panel shows de- 
rived E{B — V) versus DM1.5, plotting first the largest sample ob- 
tained for the highest position of the cut line (A(r' — Ha) = 0.03, 
in green), and then overplotting in blue the smaller sample ob- 
tained with the cut line lowered to A(r' — Ha) = 0.02. The box, 
centred on DM = 12.5 and E{B - V) = 1.15, is the result for 
NGC 7510 due to Sagar & Griffiths (1991). 



Figure 8. The 6x6 arcmin^ field centred on the test cluster 
NGC 7790, as extracted from the IPHAS database. This is the 
30-sec r' exposures from pointing 7626o. N is up and E to the 
left. 



1.3 and DM = 12.5 ± 0.3. The extinction properties of the 
candidate A dwarfs reveal that there are objects - within the 
cluster as far as we can tell - at markedly higher extinctions 
{E{B ~ V) > 1.5) than found in prior optical studies. 

The outlier grouping at DM1.5 — 11.7 that becomes 
more evident as the cut line is raised, appear to be as red- 
dened as most of the main cluster group. This similarity 
favours viewing these outliers as intrinsically brighter ob- 
jects of similar optical colour within the cluster, rather than 
as members of an independent foreground cluster. Indeed 
Piskunov et al. (2004) noticed that there was a case for a 
younger ~2 Myr-old population in NGC 7510 in addition 
to the dominant 6-10 Myr old component: the A dwarfs in 
such a component would be nearly 1 magnitude brighter as 
is implied here (e.g. the M(r') difference between 3 Myr old 
and 7 Myr old A stars is given in Table 2 as 0.9 mags). So we 
see that this selection technique begins to pick up on both 
generations of stars in NGC 7510. 

The middle panel of Fig. 7, corrected to trace A-star 
space density seemingly brings out this bimodal character 
even more. However this more physical presentation of the 
data comes at a price: the correction to space density greatly 
amplifies both the count and the Poissonian error at lower 
DM with respect to those at higher. For example for the 
data for A(r' — Ha) — 0.03 (coloured green), the error on 
the density peak in the range 11.5 < DM1.5 < 11.75 is 
nearly 40 % (just 7 stars) , while it is under 20% in the range 
12.25 < DM1.5 < 12.75 (31 stars). This needs to be borne 
in mind in the later discussion of results for Cyg OB2. 



hints that sample contamination may indeed be greater, but 
their impact is small. Another property that does not change 
very much, as the cut line is raised, is the FWHM of the 
cluster peak: it remains close to ADM = 0.4 in all 4 cases 
explored. 

As well as being clearly defined, the cluster is also at 
a mean DM1.5 and typical reddening that is in very good 
agreement with the results of Sagar & Griffiths (1991). The 
box superimposed on the reddening plot (bottom panel of 
Fig. 7) encloses the reddening range they found, and the 
error bound on the distance modulus: 1.1 < E{B — V) < 



4.3 NGC 7790 

The open cluster, NGC 7790, is beUeved to be about 100- 
120 Myrs old (Bragagha & Tosi 2006; Gupta et al. 2000) 
and hence close to the age at which A stars would begin 
leaving the main sequence. We have extracted the photom- 
etry from a 6 X 6 arcmin'^ field centred at RA 23 58.4, Dec 
-1-61 13 (J2000, see Fig. 8) and have selected candidate A 
stars following exactly the same method as for NGC 7510. 
Given the age of this cluster there is some room to question 
whether a more firmly main-sequence, slightly fainter abso- 
lute magnitude, M(r') ~ 1.9, would work better than the 
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Figure 9. Result for the open cluster NGC 7790. Derived 
E{B — V) and Z3 A/1.5, are plotted against each other. The larger 
sample shown in a lighter cyan colour was obtained from the 
31 X 31 arcmin^ overlap of IPHAS fields 7626/o, including NGC 
7790. The data shown in blue were obtained from the smaller 
6x6 arcmin^ field centred on the cluster. The cut line in both 
cases was set at A(r"' — Ha) = 0.025. The box, centred on 
DM = 12.59 and E{B - V) = 0.51, is the result for this cluster 
due to Gupta et al. (2000). 



value of 1.5 more typical of the general field. Staying with 
M(r') — 1.5, we arrive at the result shown in Fig. 9 as the 
blue datapoints. The comparison is made in this case with 
the results of the photometric study by Gupta et al. (2000) : 
they obtained E{B -V)= 0.51 ± 0.03, DM = 12.59 ± 0.15. 

Once again, it is clear this method of A-star selection 
has found a group of objects belonging to the target clus- 
ter at a range of reddenings compatible with the results 
of detailed photometric study. Had we adopted the fainter 
main sequence absolute magnitude the distribution of stars 
in Fig. 9 would still have overlapped the Gupta et al dis- 
tance, but most of the A stars would have been lying at 
lower DAI. That the brighter field absolute magnitude does 
a little better suggests that some stars have begun to evolve 
away from the main sequence. 

A second extraction across a larger field of 31 arcmin^, 
more than doubles the sample size and mildly smears the 
appearance of the E{B ~ V) versus DM1.5 diagram (the 
cyan data points in Fig. 9. But contamination of the sample 
by the general field appears slight. This outcome may imply 
that A-star selection can become a useful tool for identifying 
new clusters with economy of effort. This is the main way 
in which we use it below in the discussion of Cyg OB2. 



5 CYG OB2 

5.1 Early-A stars around the centre of Cyg OB2 

We now come on to the second main aim of this study: to 
access a lower stellar mass range than has hitherto been 
open to examination in Cyg OB2 and its environment. The 
reward for doing this will be, at minimum, a better picture 
of the angular extent of this association. 
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Figure 10. The sampled variation of pA versus provisional dis- 
tance modulus (DA/i = dereddened r' magnitudes less an 
adopted median absolute magnitude of 1.0) for fields 6010 and 
5985 to either side of the centre of Cyg OB2. In each panel the 
four curves shown are for A(r' — Ha) = 0.015, 0.02, 0.025 and 
0.03, and are again respectively black, blue, red and green, as in 
Fig. 7. The adopted bin size is A{DM) = 0.25. 



We first consider the selection and dereddening of can- 
didate early-A dwarfs in IPHAS fields 6010 and 5985, that 
straddle the centre of Cyg OB2 (Fig. 1). The colour-colour 
plane for field overlap 6010 was shown as Fig. 2, from which 
it was noted that Cyg 0B2 is apparent as an excess of ob- 
jects in the approximate reddening range 4.5 < Ay < 7. As 
in Fig. 2, no selected object will be fainter then r' = 20. 
As for NGC 7510 and NGC 7790, this cut-off magnitude is 
2 magnitudes brighter than the reduction pipeline estimate 
of the magnitude limit of the observations. In this and all 
other respects, the method of selection is the same as those 
applied to these test-case clusters. 

In Fig. 10, we present the space density distributions for 
provisional distance moduli, DM\.o, obtained from early-A 
dwarf candidates for different offsets of the IPHAS limit line 
for both fields: as the offset is increased, so the sample size 
increases. We choose M{r') — 1.0 for now, noting that this 
corresponds to an age in the region of 5 Myr. In its lowest 
position the cut line is placed at A{r' — Ha) = 0.015 above 
the AOV reddening line. The increasing contamination with 
increasing offset potentially brings in intrinsically redder, 
and fainter, stars, that might bias the derived distribution 
toward over-estimated reddenings and underestimated dis- 
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Figure 11. Visual extinctions versus provisional distance moduli 
DMi.o for fields 6010 and 5985 straddling the centre of Cyg OB2. 
The dark blue points correspond to objects selected below the 
A(r' — Ha) = 0.02 cut line (cf Fig. 10), while the lighter green 
points are the objects added to the sample for the highest cut line 
(A(r' — Ha) = 0.03). The box drawn in black and the green line 
are discussed in the text. The box in blue delineates a range in 
DM\x) also picked out in figs 13 and 14. 



tance moduli. A very small effect, possibly due to this, may 
be apparent in the populous, less reddened field (5985) to 
the NW of the nominal Cyg OB2 centre - while the over- 
all sample size grows from 142 objects (black curve) to 223 
(green curve). In field 6010 the shift is negligible and per- 
haps just positive, while the sample size increases from 37 
objects (black curve) to 77 (green curve). To limit Poisso- 
nian noise, we will prefer the largest samples obtained on 
setting A(r' - Ha) = 0.03. 

The striking feature of Fig. 10 is the large difference 
in the A-star density distributions contained in the two 
neighbouring fields. For 6010, the peak in the range 11.0 < 
DMi.o < 11.25 sits at 3-sigma contrast with respect to the 
less well-populated bins to either side of it. No such peak 
is apparent in 5985 and the overall population seems less 
dense and more distant. To better understand the origin of 
this, we show in Fig. 11, the relationship between reddening 
(Av) and DMi.o in both cases. In field 6010 just to the SE of 
the nominal centre of Cyg 0B2, it is evident that the stars 
making up the peak population around DMi.o — H are 
seen through visual extinctions of 5 to 8 magnitudes. This 
contrasts with 5985 where the most common reddening to 



DMi.o ~ 11 is 4 to 6 magnitudes. The other clear point 
of contrast is that it is the slower rise of reddening with 
distance within 5985/o that allows us sight, in this field, 
of so very many candidate A dwarfs at DMi.o > 12. We 
are seeing these objects through a relative reddening hole 
that has been known about since the work of Reddish et al. 
(1967, see also Dickel & Wendker 1978, Knodlseder 2000). 
It is also evident (from Fig, 11), that the r' = 20 magnitude 
limit applied slices through this more distant population. 

The box drawn in black in both panels of Fig. 11 spans 
the range in distance modulus and extinction obtained by 
Hanson (2003) in her study of brighter OB stars. The mean 
DAI she favoured was 10.8. The vertical green line in both 
panels marks DM = 11.2, obtained by Massey & Thomp- 
son (1991). The tidy picture seen in the case of NGC 7510, a 
relatively isolated young cluster, is not repeated here. In the 
upper panel, for 6010, the Hanson (2003) box looks to be dis- 
placed to lower distance moduli than the selected A-dwarfs: 
to bring the two into alignment it would be necessary to 
increase the assumed A-dwarf absolute magnitude by ~0.7. 
The Massey & Thompson (1991) distance fares rather better 
in that the offset is only about 0.3 magnitudes. With either 
of these comparisons there is a problem in that increasing 
the adopted absolute magnitude amounts to raising the age 
of the population, when existing estimates for the age of 
Cyg OB2 are less than 5 Myr. 

In 5985 (lower panel of Fig. 11), the Hanson (2003) 
Cyg OB2 distance and reddening range exhibits even less 
connection with the apparent distribution of candidate 
early-A stars. Much of the early-A population in this field is 
unlikely to belong to any cluster at a distance below ^^2 kpc 
{DMi.o 11.5). Indeed, there may be a steady rise in red- 
dening with increasing distance modulus in this field, with 
little imprint of Cyg OB2 itself. 

5.2 Candidate A dwarfs across the wider Cyg 
OB 2 region 

The larger picture in this region is shown in Fig. 12. Dis- 
tributions are only shown for fields in which there are more 
than 50 candidate early-A dwarfs (for A(r' — Ha) — 0.03). 
The number of candidate objects per field was given in Ta- 
ble 1. The better populated fields lie mainly in the cen- 
tre and west of the region studied. The outstanding feature 
of the A-star density profiles as a function of DA/i.o, ap- 
parent from Fig. 12, is the repeated density peak in the 
11.0 < DMi.o < 11.25 bin in the southern fields (6010, 
5976, 6035, 6001 and 5985). High densities at similar dis- 
tance moduli are also apparent in 5951 and 5998 - but not 
in the two fields immediately to the north of the nominal 
Cyg OB2 centre, 6021 and 5985. 

To show in greater detail what this may mean, we plot 
reddenings against DA/i.o for fields 6001 and 5976 in Fig. 13. 
The blue box is drawn around the main concentration of can- 
didate early-A dwarfs found here (10.8 < DMi.o < 11-6). 
Objects in this range make up the overwhelmingly domi- 
nant pA peak seen in Fig. 12 for field 6001 (39 stars, but 
note that the pA peak closer to DMi.o ~ 10 is created by 
just 5 stars and so is much less significant). The reddening 
in these fields is somewhat higher than in the central part of 
Cyg OB2 (6010) to the north: for most stars here, Av > 6. In 
field 5976 there is a tail of detected objects to much greater 
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Figure 12. Derived candidate early-A dwarf densities as a func- 
tion of DAfi.o for fields in and around Cyg OB2. Tlie relative 
positioning of the histograms roughly mimics their on-sky posi- 
tioning, as shown in Fig. 1. Hence, roughly speaking, N is up 
and E to the left. To increase the sample sizes, the cut line in 
the IPHAS (r' — Ha, r' — i') plane has been set at maximum 
(A(r' — Ha = 0.03). The most nearly central field, 6010, is picked 
out in blue. Overlap fields listed in Table 1 with fewer than 50 
candidate early-A dwarfs are not shown. In every panel, the axis 
ranges are the same and, in some, a vertical line is drawn at 
DMi.o = 11.1. 
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Figure 13. Visual extinctions versus provisional distance moduli 
DA'Ii o for fields 6001 and 5976 to the south-west of the centre 
of Cyg OB2. The dark blue points correspond to objects selected 
below the A (r' - Ha) = 0.02 cut line (cf Fig. 10) , while the lighter 
green points are the objects added to the sample for the highest 
cut line {A(r' -Ha) = 0.03). The DMi.o range picked out by the 
same blue box in each panel is the peak density region apparent 
from Fig. 12. To ease comparison, this region is also identified in 
both Figs. 11 and 14. 



incompleteness of this kind is much more of an issue in the 
extreme north of the region. 



distances (DMi.o > 12) which is absent from 6001, suggest- 
ing that the typical extinction affecting sightlines in 5976 is 
less severe. This also raises the possibility that the greater 
obscuration in field 6001 induces more sample incomplete- 
ness: undercounting of A stars in this field at DMi.o H-O 
may be more significant. 

Finally, we take a look at 0.5-1 degree N and NW of the 
centre of Cyg OB2, at fields 5998 and 5962 (Fig. 14). Field 
5998/o is rather similar to 6010, and fields further south. 
The character of 5962 is quite difi'erent: there is not a con- 
centration of stars in the DMi.o range, 11.0 - 11.5, and the 
overall number of candidate A dwarfs is relatively low. The 
modest pile up at Av ~ 8 at distance moduli 12.0 to 13.0 
in 5962 (Fig. 14) suggests there may be more undiscovered 
at still higher reddenings, fainter than the imposed r' = 20 
limit. Inspection of similar reddening plots (not shown) for 
the fields neighbouring 5962 supports the impression that 



5.3 Implications for the extent and nature of 
Cyg OB2 

The main implication of Fig. 12 is the clear A-star density 
peak centred on 11.0 < DMi.o < 11.25 in the southern 
fields, that is most pronounced in field 6001. Is this the real 
Cyg 0B2? This can be looked at in a different way by map- 
ping the positions of the stars with dereddened r' magni- 
tudes falling in the range 11.75 < v'q < 12.5 that contribute 
to this apparent 'clustering' (Fig. 15). 

On the basis of 2MASS data, Knodlseder (2000) demon- 
strated a peak in stellar density at NIR wavelengths that 
coincides with the heart of Cyg 0B2 itself. This region is 
enclosed by the the red box drawn in Fig. 15. There is not 
a matching peak in the concentration of A stars here: the 
IPHAS field pair 6010 does not stand out in this way. More- 
over, there is a relatively empty region cutting across the NE 
corner of 6010, continuing through the south of adjacent 
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Figure 15. The spatial distribution of candidate A dwarfs in the vicinity of Cyg OB2, with dereddened r' magnitudes in the range 
11.75 < < 12.5. The area bounded by the red line roughly matches the area of peak NIR, stellar density as identified by Knodlseder 
(2000). The area bound by the blue line picks out the main concentration of optically identified OB stars (Massey & Thompson 1991). 
In the NE the green box encloses the stellar population associated with DR 21 discussed by Kumar et al. (2007). 



fields 6021 and 5985. This underpopulated region is mucli 
tiie same as that occupied by the OB stars picked out in 
Massey & Thompson's (1991) optical study (the rough out- 
line to this region is picked out in blue in Fig. 15, see also 
Fig. 4 of Hanson 2003). Inspection of the original IPHAS 
imaging confirms that this anti-correlation is not due to the 
OB stars creating a bright background swamping the fainter 
A stars. 

As shown in Fig. 15 the east of the region is even more 
striking for the lack of candidate A dwarfs. This absence was 
already foretold by the data in Table 1, showing that there 
are few ordinary A stars to be had at any distance mod- 
ulus. The appearance of a few in the extreme north-east, 
in the vicinity of DR 21 (marked in green), may signal the 
emergence of the new generation of stars linked to this well- 
known HIT region. The median rg for the group is 12.33. If 
indeed they are as young (1-3 Myrs) as the population of ob- 
jects analysed by Kumar et al.(2007), they would be 2.9 kpc 
away (adopting M{r') = 0) - the distance also favoured by 
Kumar et al. (2007). 

The major concentration of A dwarf candidates in 
Fig. 15 is in the southern fields that also stood out in Fig. 12, 
albeit with a zone of avoidance in the south west (lower part 
of 5976, extending to most of 5939 further west). If the cen- 
tral field, 6010, is included with these fields (5951, 5965, 



5976, 6001, 6035), the total number of plotted objects is 
239 - an average of 40 per field. 6001 is the most populous 
with 49 stars, 6010 the least with 32. This leaves just 125 
objects to be shared among the remaining 14 fields in this 
magnitude range. In broad brush terms, we have roughly 
200 A stars spread across an area about 1 degree across, 
distributed around RA 20h 32m and Dec +41 00 (2000). 

The typical dereddened magnitude for this grouping is 
r'o ~ 12.1. To be sure of the distance to this structure, an 
age must be known. The assumption made for Fig. 12 was 
consistent with an age of 5 Myr, and placed this clustering at 
distance of 1.7 kpc (DMi.o — 11.1). If the typical age can be 
allowed to be as low as 2 Myr, then M(r') is required to be 
as bright as 0.0, driving the distance up to an improbable 
2.6 kpc (DA/o.o = 12.1). A further argument against this 
solution is that we are seeing evidence here that such young 
A star populations tend to be sparse - as illustrated by 
the low counts linked to both DR 21 in the NE and to the 
main central concentration of optically-detected OB stars - 
both likely to be only 1-3 Myrs old. It is more plausible 
to consider ages greater than 5 Myrs, rather than less. For 
example, an age of ~7 Myrs imposes M(r') — 1.4 (Table 2) 
which then points to a distance modulus close to the value of 
10.8, equivalent to a distance of 1.4 kpc, that was favoured 
by Hanson 2003 based on revised OB-star calibrations. 
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Figure 14. Visual extinctions versus provisional distance moduli 
DMi.o for fields 5998 and 5962 to the north and north-west of 
the centre of Cyg OB2. The colours have the same significance as 
in Figs. 11 and 13. 

Accepting an age spread in and around Cyg OB2 seems 
the best option. If there are very young A-type objects em- 
bedded among the OB stars, they are few in number and/or 
escaping selection by virtue of either modest Ha emission or 
high extinction {Av > 10). The larger numbers of A stars 
to the south of Cyg 0B2 centre give away the presence of 
the larger, somewhat older, main cluster. In effect Cyg OB2 
can be seen as a further example of cluster 'rejuvenation', 
but on a much grander scale than that already noted in 
the test case of NGC 7510. What is unclear at the moment 
is whether the rejuvenation is on the NE periphery of the 
main cluster, or whether it is more centrally located within 
a cluster that extends as far north as field 5998, which also 
exhibits an A-star density peak in the appropriate dered- 
dened magnitude range (Fig. 12). The less elegant solution 
at this point is to view the swathe of A stars to the south 
of the nominal centre of Cyg 0B2 as evidence of a distinct 
cluster, of unknown age and distance, unconnected with the 
known OB stars. 

If we can associate ~200 A0-A5 near main sequence 
stars with a 5-7 Myr old southerly cluster, then the com- 
bination of Siess et al. (2000) isochrones for this age range 
and the Kroupa (2001) flattening of the Salpeter IMF be- 
low 0.5 M0 would imply a total cluster mass in the region 
of 10,000 Mq, and an OB star count of around 200. There 
are a number of respects in which this estimate is surely on 



the low side. First, a minor point is that the sky coverage 
is only 90 percent. But more important is the fact that the 
selection of A stars applied here has been fairly conservative 
in order to keep sample contamination minimal: had the cut 
line been raised to A(r' — Ha) = 0.05, the number of can- 
didate A dwarfs in the region would nearly double. There is 
also clear evidence of sample incompleteness due to high and 
patchy reddening - a problem that appears more marked in 
the more southerly fields 6001 and 5976 with higher pA than 
near the OB star concentration. And finally, this accounting 
does not include the immediate environs of the OB stars in 
which fewer stars of early-A colour have yet appeared: ap- 
plying the same approach as above for a star count of ~30 
(for e.g. field 6010) at the much younger age of ~2 Myrs 
adds a further ~10,000 M0. A final proper accounting of 
the total mass in the region of 30,000 to 40,000 Mq would 
not be surprising. 

The reputation of Cyg OB2 as the outstanding visible 
northern OB association remains untarnished. However, the 
upper limit of the claim by Knodlseder (2000) is difficult to 
sustain: if this region truly is home to over 2000 OB stars 
then this appraisal of the A star content has under-counted 
by an order of magnitude. At least in the core region of Cyg 
OB2 where the reddening is less severe, so large a discrep- 
ancy is unlikely. 



6 CLOSING REMARKS 

This study has had two goals. The first of these was to show 
how the IPHAS {r' — Ha, r' — i') plane may be exploited 
in a very simple way to yield reliable samples of near main- 
sequence A0-A5 stars. Such objects in the general field, se- 
lected as described here and validated by HIPPARCOS data, 
will present themselves with a typical absolute magnitude 
M(r') of 1.5 - 1.6. Fortunately, for stellar ages of 10 Myrs 
or more there is little change in this absolute magnitude. 
For simple single-aged populations a magnitude spread of 
FWHM ~ 0.5 is to be expected. Demonstrations of this 
selection at work have been provided by the examples of the 
open clusters NGC 7510 and NGC 7790. 

The second aim has been to apply this A-star selection 
technique to the challenging example of Cyg OB2. This spec- 
tacular northern cluster, immersed in the Cygnus X region 
of the Galactic Plane, continues to be debated in the liter- 
ature on account of its uncertain dimensions and distance. 
The new capability directly bestowed by IPHAS photometry 
(unsupported, thus far, by spectroscopy) is that of identify- 
ing the on-sky distribution of and reddenings of over 1000 
A0-A5 stars towards, around, within, and beyond Cyg OB2. 
The result of this is a clear peak in A-star density, compris- 
ing ~200 objects, roughly centred on RA 20h 32m and Dec 
+41 00 (2000), some 20 arcmin south of the VI Cyg No.8 
trapezium normally taken as the centre of Cyg 0B2. 

The most parsimonious interpretation of the results ob- 
tained is to treat the southern clustering of A stars and the 
already known OB star concentration as part of the same 
cluster, compatible with a distance of 1.4 kpc {DM ~ 10.8, 
Hanson 2003) if the A stars are ~7 Myr old, or with 1.7 kpc 
{DM = 11.2, Massey & Thompson 1991) if their typical 
age is ^^5 Myr. Either solution requires some age contrast 
between the A stars (5-7 Myrs old) and the main group of 
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OB stars (1-3 Myrs old, Hanson 2003). This may amount 
to two main bursts ol star formation, or to an indication 
of a substantial age spread (cf. the case of the ONC clus- 
ter discussed by Palla et al. 2005). The presence of a range 
of ages may help to explain why Cyg OB2 is not associ- 
ated with bright H ll emission despite its very evident rich 
O star population. The massive cluster, R 136 in the Large 
Magellanic Cloud exhibits similar behaviour in that its inter- 
mediate mass stars appear to be ~6 Myrs old, as compared 
with 1-2 Myrs for the most massive stars (sec Massoy 2003) . 
Indeed, in Cyg OB2 itself, signs of non-coeval star forma- 
tion have been noticed before (Massey & Thompson 1991, 
Hanson 2003). The mass linked to the newly- found A stars 
in the cluster is at least 10,000 Mq. As a cautionary note, 
we point out that it remains possible (if less attractive, see 
Schneider et al 2006) that the main A-star clustering is at a 
different distance from the OB stars. 

There are two wider applications of this new technique 
for selecting A stars from IPHAS. First, it is a new tool for 
use in exploring the structure of the Galactic thin disk. This 
will be especially incisive when applied to the less reddened 
outer Galaxy. We have not discussed metallicity here, but it 
should not create much difficulty: on the one hand, the weak 
line-blanketing at red wavelengths implies that the r' — i' 
colour of A stars will not be metallicity-sensitive, and, on the 
other, the evidence to date is that metallicities {[Fe/H]) in 
the thin disk at large Galactic radii are not greatly subsolar 
(e.g. Carraro et al 2007, who find a mean [Fe/H] of ^ —0.35 
for Rg between 12 and 21 kpc). Second, A-star selection over 
large sky areas allows a broad census to be conducted of the 
early evolution of A stars by, for example, identifying the 
fraction with clear infrared excesses. A study of this kind 
has already been undertaken by Hales et al. (2008), who 
have combined IPHAS photometry with Spitzer/GLIMPSE 
and 2MASS data. 

A work for the future is to generalise the A-star selec- 
tion scheme to begin to extract a broader range of spectral 
types from the main stellar locus in the IPHAS colour-colour 
plane. In principle it is possible to map the positions of A- 
K near main-sequence stars as a function of reddening onto 
unique (r' — Ha, r — i') values - if an effective means to dis- 
tnguish and separate out higher luminosity classes can be 
identified. An exploration of this opportunity is now under- 
way (Sale et al., in preparation). 
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